Interaction of vortices in thin superconducting films and 
Berezinskii-Kosterlitz-Thouless transition 
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The precondition for the BKT transition in thin superconducting films, the logarithmic intervortex 
interaction, is satisfied at distances short relative to A = 2A^/d, A is the London penetration depth 
of the bulk material and d is the film thickness. For this reason, the search for the transition has 
been conducted in samples of the size L < A. It is argued below that film edges turn the interaction 
into near exponential (short-range) thus making the BKT transition impossible. If however the 
substrate is superconducting and separated from the film by an insulated layer, the logarithmic 
intervortex interaction is recovered and the BKT transition should be observable. 



PACS numbers: 74.76.-w 

Introduction. The seminal prediction by Berezinskii 
[H, Kosterlitz, and Thouless (BKT) of the phase tran- 
sition caused by fluctuations-induced topological exci- 
tations in two-dimensional (2D) systems has been con- 
firmed in a number of experiments with superfluid films. 
Many attempts to do the same with thin superconducting 
films were less convincing, for a review of early theoreti- 
cal and experimental work see, e.g., Ref-i Later trans- 
port data on superconducting films remainedj.nconclu- 
sive 0, @| ■ The recent transport measurements on high 
quality ultra-thin high-T^ superconducting films failed to 
show the jump in the exponent a of current- voltage char- 
acteristics, V (X I'^, a signature of the BKT transition 
Q . The discussion persists up to this day [^, Q . 

The precondition for the BKT transition, the logarith- 
mic intervortex interaction, is satisfied only at distances 
short relative to A = 2\^/d (A is the London penetration 
depth of the bulk material and d is the film thickness). 
For this reason, the search for the transition has been 
conducted in samples of the size L < A. It is argued here 
that even in small samples the interaction is not logarith- 
mic: due to the boundary conditions at the film edges the 
interaction turns into a short-range near exponential de- 
cay. The main conclusion therefore is that the BKT tran- 
sition could not happen in thin superconducting films of 
any size on insulating substrates (unlike in layered com- 
pounds where the interaction between pancake vortices 
is logarithmic tlDl])- If however the film is placed on a 
superconducting substrate and separated from it by an 
insulating layer of the thickness s, the logarithmic inter- 
action on all distances greater than s is recovered and 
the BKT transition should take place. 

The situation is different for ac response of thin 
films. The characteristic separation of vortices 1^^ cxl/oj 
contributing to the 2D response might be small at large 
frequencies lu and therefore can be recorded even in 



small samples, see Ref.[ll| and references therein. We do 
not consider ac phenomena here. 

Approach. We begin with a brief review of an approach 



to vortices in thin films suggested in Ref. 12]; although 
not common this approach allows one to evaluate ener- 
gies in a direct manner, the adva ntag e relevant for our 
purpose. As was stressed by Pearl [13[, a large contribu- 
tion to the vortex energy in thin films comes from stray 
fields. In fact, the problem of a vortex in a thin film is 
reduced to that of the field distribution in free space sub- 
ject to certain boundary conditions at the film surface. 
Since curl h = div h = out the film, one can introduce 
a scalar potential for the outside field: 



(1) 



To formulate the boundary conditions for the outside 
Laplace problem, consider a film of thickness d <C A oc- 
cupying the xy plane. For a vortex at r = 0, the London 
equations for the film interior read after averaging over 
the film thickness: 



27rA , ^ ^, , 
hz H curing = 0o 5{r) , 



(2) 



where g{r) is the sheet current density, r — {x,y) and 
(/)o is the flux quantum. 

In a thin film, the Maxwell equation Anj/c = curl/i is 
reduced to relations between the sheet current and dis- 
continuities of the tangential fields: 

47rt//c = £ X [ h{+0) - h{-0)] . (3) 

One substituts Eq. ([3]) in ((21) and uses div h ~ Oto obtain: 

hz ~ A[dzhz{+0) - dzhz{~0)]/2 = 0o'5(r) . (4) 

This equation expressed in terms of the potential ip along 
with conditions at infinity constitute the boundary con- 
ditions for the Laplace problem, Eq. ([T]), for the field dis- 
tribution outside the film. 

Consider the case in which the half-spaces above and 
under the film are vacuum (or a non-magnetic insulator) . 
The general form of the potential that vanishes at z ^ 
+00 of the empty upper half-space is 
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Here, (p{k) is the 2D Fourier transform and fc^ — k'^ + k^. 
In the lower half-space one has to replace z ^ —z. Due 
to the symmetry, the boundary condition ^ becomes 
d^ip{+0) - Ad'^ip{+0) = (poSir) that yields the Fourier 
transform of the potential in the upper half-space: 



(p = (fi + (p2 and 



(9) 



^ = 



fc(l + fcA) 



(6) 



This provides hz{k) = —kip{k) and hx^y{k) — ik^^y ip{k), 
i.e. the fields everywhere in the free space as well as the 
sheet currents due to a single vortex in an infinite film. 
One can easily verify that this solution coincides with 
that given by Pearl [l3|. 

Energy. We now establish connection between vortex 
energy and the potential ip. We begin with the general 
situation of a vortex in a finite bulk sample. The energy 
consists of the London energy (magnetic -I- kinetic) in- 
side the sample, e^*) = J[h^ + {An Xj / c)^]dV / Sn and the 
magnetic energy outside, e^"^ = / h'^dV/Sir. Then, for 
the potential introduced in Eq. ^ and gauged to zero 
at CO (which is possible in zero applied field) one has 
87re*^°^ = / (fih ■ dS where the integral is over the sample 
surface with dS directed inward the material. The Lon- 
don part is transformed integrating by parts the kinetic 
term: Sttc^'^ — (47rA^/c) ^(j x where the integral is 
over the sample surface and the surface of the vortex core. 
The integral over the sample surface is further trans- 
formed: § dS-ijxVtp) = f dS-ip{V xj) [l3|. Combining 
the result with e^°'\ one obtains <f dS-ip{h + 4:nX'^cnTlj/c) 
where the expression in parenthesis is <l)ov6'^^^{r — r„) 
where v is the direction of the vortex crossing the sur- 
face at the point r„, and S^'^^ (r — ry) is the 2D S function. 
One then obtains: 



in Ac 



dS-{hxj); (7) 



Tent and Vex a^G the positions of the vortex entry and exit 
at the sample surface (the vortex is assumed to cross the 
sample surface at right angles; otherwise, one should mul- 
tiply the potentials by cosines of corresponding angles). 
For more than one vortex one has to sum up expressions 
([7]) over all vortices. If more than one superconductor is 
present, Eq. ^ still holds, but </? as a solution of Laplace 
equation outside superconductors will be affected by the 
presence of all superconductors. 

For thin films, the integral over the core surface (cx d) 
can be neglected: 



00 



(8) 



where the notation Dipi, = (pi,{rent) — ^vi^ex) for the 
I/- vortex is introduced for brevity. Due to linearity of the 
Laplace and London equations, one has for two vortices 



where the arguments 1,2 are positions of vortices. 
Clearly, the self-energy of the first vortex is 



and the interaction energy is given by 



Stt 



(10) 



(11) 



Infinite film in vacuum. Equations (jlO[) and (|11[) are 
quite general and hold for films of any lateral size. In 
particular, for an infinite film in vacuum one obtains with 
the help of the potential ([5]): 



.(0) 



47r 



d^k 



(27r)2 k{l + kK) 



(12) 



where the cutoff at fcmax ~ 1/^ is introduced to a loga- 
rithmically divergent integral. Similarly, one finds: 



SttA 



Ho 



A 



Yo 



A 



(13) 



where Hq and Yq are Struve and Bessel functions and 
r is the distance between vortices. At distances r ^ A 
this yields Cmt = '/'o/^"'^^ (as for two point "magnetic 
charges" ^o/Stt) showing that at large distances the role 
of stray fields is dominant. For r <C A, the interaction is 
logarithmic. 



Att'^A r 



(14) 



that led to a statement that the BKT transition should 
be searched for in samples of the size L <^ A. We show 
below that this statement is incorrect. 

Small samples in vacuum. To consider small samples, 
one turns back to the basic Eq. The currents g{r) can 
be found by solving Eq. ^ combined with the continuity 
equation and the Biot-Savart integral which relates the 
field h-, to the surface current: 



divg = 0, K{r) = j d^r'[g{r') x R/cR^ 



(15) 



R = r — r' . To satisfy divg = it is convenient to deal 
with a scalar stream function G{r) such that g — curlGi. 
Alternatively, the sheet current can be expressed in terms 
of the derivatives of the potential <f{r, +0) at the upper 
film face: 



(16) 
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Therefore, G{r) is proportional to Lp{r, z = +0): 
ip{r,+0) = -27rG(r)/c 



(17) 



(a possible additive constant is zero since both cp and 
G are gauged to zero at oo). Since the energies are ex- 
pressed in terms of (p at vortex positions, Eq. p7|) shows 
that to evaluate vortex energies and their interaction it 
suffices to find the current distribution G(r). 

The problem of current distribution for samples of arbi- 
trary size is difficult: the London Eq. ^ combined with 
Biot-Savart law (jlSp make an integro-differential equa- 
tion for g. However, for small samples, L <g; A, it is man- 
ageable because as is seen from Eq. ([T5|) the field ^ g/c 
whereas the term with current derivatives in Eq. ^ is of 
the order Ag/cL. Then, Eq. ([2]) can be rewritten as a 
Poisson equation for G: 



V^G -(c(/)o/27rA)(5(r - a) 



(18) 



The normal to the edge component of the current 
vanishes, i.e., G is a constant along the edges which 
can be set zero. Thus, G is equivalent to the elec- 
trostatic potential of a linear charge q = c^io/Sir^A 
situated at the vortex position a and parallel to the 
side surface of the grounded metal cylinder with the 
crossection coinciding with the thin-film sample. A 
rich library of the 2D electrostatics is instrumental 
in solving for vortex currents for a variety of film 
shapes with linear dimensions less than A; see, e.g., 
Ref. [3l for the solution of Eq. (fT8|) for a rectangular film. 



Thin-film strips in vacuum. Most of experiments of 
our interest are done on long thin-film strips. For a strip 
along y with edges at a; = and x ^ W A, the solution 
of Eq. (|18p for a vortex placed at a = [qx, ay) with zero 
boundary conditions at the edges reads: 



tanh 



sin(7raa;) sin(7ra;) 



G 

2q cosh[7r(2/ — ay)] ~ cos(7raa;) cos(7ra;) ' 



(19) 



where for brevity W is used as a unit length (see, e.g., 
Refs. 15[) or [16|. Due to the edges, G does not de- 
pend exclusively on |r — a|; note however that G(r, a) = 
G(a, r). It is straightforward now to find the energies of 
vortices in a strip and their interaction with the help of 
general formulas (HO]), (HIl), and ([TTll . 

The self-energy of a vortex at the position x = a, y = Q, 
e = (/>oG(a)/2c, is logarithmically divergent. Introducing 
the cutoff at |r — a| — ^ one obtains [12|: 



(20) 



The interaction (|lip of two vortices at Oi and 02 is 
(1)1 cosh[7r(yi - ya)] - cos[7r(a;i + X2)] 



In 

87r2A cosh[7r(yi - 7/2)] - cos[7r(a;i - X2)] 



(21) 



where the coordinates are given again in units of W. 
(Formally similar interaction exists between Abrikosov 
vortices parallel to the thin film [l3|-) Hence the interac- 
tion is not proportional to the logarithm of the intervor- 
tex distance as required by BKT (except at distances of 
the order ^ where the present theory breaks down). E.g., 
for yi = 2/2 we have (in conventional units): 



47r2A 



In 



sin[7r(a;i + X2)/2W] 



sin[7r(a;i - X2)/2W] 



(22) 



For vortices in the strip middle {xi = X2 = W/2) sepa- 
rated by y, 



I- In 



Att^A 

This gives for y > W: 



coth 



27r2A 



exp 



\2WJ 



27ry\ 
W J 



(23) 



(24) 



Thus, for separations larger than W/2tt, the intervortex 
interaction is in fact exponentially weak and short-range. 
Physically, this happens because the second vortex feels 
not only the "bare" first, but also the infinite chain 
of ± images which one has to introduce to satisfy the 
boundary conditions at the edges. The possibility of 
BKT transition in thin-film bridges of YBCO similar to 
those studied in has been examined experimentally 
3 and theoretically 0] (see references therein) with no 
positive outcome. Given the above argument, one can 
say that the short range interaction in small samples is 
the reason for this failure. 

Superconducting substrate. Up to this point, we have 
discussed a "free standing" film or, better to say, a 
thin superconducting film placed on an insulating non- 
magnetic substrate so that the symmetry of the upper 
and lower half-spaces could be utilised. The situation 
changes if the film is placed on a superconducting sub- 
strate being separated from it by an insulating layer to 
prevent the Josephson coupling to the substrate. The 
vortex magnetic flux (j>Q is channeled into the space be- 
tween the film and the screening substrate so that the 
radial field component along with azimuthal sheet cur- 
rents vary as 1/r at large distances. This results in the 
Lorentz interaction force of two vortices varying as 1 /r 
and, therefore, in the logarithmic interaction energy. 

To confirm this quantitatively, consider a film situated 
at a z = s above a superconducting substrate occupying 
the half-space z < 0. If s is large enough to suppress the 
Josephson interaction, the problem is reduced to solving 
the Laplace equation ([T]) for the potential fAix, z) in the 
free space z > s (the domain A) and for ipsif, z) in the 
domain B between the film and the substrate, < z < s. 
A simple method of doing this is described in Ref. 1^. 
One looks for ipA{f, z) in the form ([5]) with the 2D Fourier 



4 



transform (pA{k) e Similarly, in the finite domain B, 
(/3b(A;,2) = Cl(fc)e'=^+C2(fc)e-'=^ The boundary condi- 
tions at z = s are given by the field continuity and by the 
London equation ^ . One further simplifies the problem 
of the substrate screening by setting the substrate pene- 
tration depth to zero, i.e., hz{z — 0) — at the substrate 
surface; this gives Ci ~ C2 = C. After simple algebra 
one obtains: 



non-superconducting substrate. On the contrary, if the 
substrate is superconducting and separated from the film 
by an insulated layer of a thicknes s, the interaction is 
logarithmic at all distances r > s and the BKT transition 
should be observable in a dc type of experiment; e.g., the 
dc power-law current-voltage characteristics should show 
the well-known jump of the power exponent at the BKT 
transition. 



2Ce'"' sinh ks 



(po e sinh ks 



fc(sinh ks + kAe'"'/2) ' 



(25) 



Given the potentials above and under the film, one uses 
Eq. pO)) to calculate the vortex energy: 



e = 



dk e 



ks 



167r2 J sinh ks + kA e'"' /2 ' 



(26) 



The lower limit of integration is the inverse sample size 
whereas the upper one is 1/^. One can estimate this 
integral splitting the integration domain in two: s/A < 
ks < 1 and 1 < ks < s/^. The contributions to the 
integral are estimated as (2/A) \n{L/s) and (4/A) ln(s/^). 
This gives 



47r2A 



'° In.^ 



(27) 



The interaction energy is found with the help of the 
general result ([TT|) : 



, , cosh ks -I- e*"'* sinh ks , , 
sinhfc.-KMe^V2 ' ^''^ 



The oscillating Bessel function truncates the integral at 
k K, \/r^ so that in the relevant part of the integration 
domain ks < s/r for intervortex distances exceeding 
s. One readily estimates: 



47r2A r 



r > s . 



(29) 



Thus, the superconducting film situated parallel to the 
surface of a bulk superconductor should exhibit the 
BKT transition, unlike the case of insulating substrate, 
a verifiable conclusion. 

To summarize, it is shown that in small thin-film sam- 
ples on insulating substrates, edge effects modify the 
vortex-vortex interaction making it short-range, unlike 
the logarithmic long-range interaction needed for the 
BKT transition. In large free-standing films this tran- 
sition cannot happen because of the l/r interaction via 
the stray fields. This makes the BKT transition impos- 
sible in thin films of any size if they are supported by a 



The idea of this work emerged while studying the IV 
characteristics of thin-film bridges taken by F. Tafuri and 
J. Kirtley. Also, the author is indebted to L. Bulacvskii 
and R. Mints for helpful discussions. The work is sup- 
ported by the Office of Basic Energy Sciences, US DOE. 
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